\(\int \frac {(a+b \tan (e+f x))^2 (A+B \tan (e+f x)+C \tan ^2(e+f x))}{(c+d \tan (e+f x))^2} \, dx\) [78]

Optimal result

Integrand size = 45, antiderivative size = 417 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=-\frac {\left (a^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-b^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-2 a b \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) x}{\left (c^2+d^2\right )^2}+\frac {\left (2 a b \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )+a^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )-b^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \log (\cos (e+f x))}{\left (c^2+d^2\right )^2 f}-\frac {(b c-a d) \left (b \left (2 c^4 C-B c^3 d+4 c^2 C d^2-3 B c d^3+2 A d^4\right )+a d^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \log (c+d \tan (e+f x))}{d^3 \left (c^2+d^2\right )^2 f}+\frac {b^2 \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan (e+f x)}{d^2 \left (c^2+d^2\right ) f}-\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))} \]

[Out]

-(a^2*(c^2*C-2*B*c*d-C*d^2-A*(c^2-d^2))-b^2*(c^2*C-2*B*c*d-C*d^2-A*(c^2-d^2))-2*a*b*(2*c*(A-C)*d-B*(c^2-d^2)))
*x/(c^2+d^2)^2+(2*a*b*(c^2*C-2*B*c*d-C*d^2-A*(c^2-d^2))+a^2*(2*c*(A-C)*d-B*(c^2-d^2))-b^2*(2*c*(A-C)*d-B*(c^2-
d^2)))*ln(cos(f*x+e))/(c^2+d^2)^2/f-(-a*d+b*c)*(b*(2*A*d^4-B*c^3*d-3*B*c*d^3+2*C*c^4+4*C*c^2*d^2)+a*d^2*(2*c*(
A-C)*d-B*(c^2-d^2)))*ln(c+d*tan(f*x+e))/d^3/(c^2+d^2)^2/f+b^2*(2*c^2*C-B*c*d+(A+C)*d^2)*tan(f*x+e)/d^2/(c^2+d^
2)/f-(A*d^2-B*c*d+C*c^2)*(a+b*tan(f*x+e))^2/d/(c^2+d^2)/f/(c+d*tan(f*x+e))

Rubi [A] (verified)

Time = 1.20 (sec) , antiderivative size = 417, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.133, Rules used = {3726, 3718, 3707, 3698, 31, 3556} \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\frac {\log (\cos (e+f x)) \left (a^2 \left (2 c d (A-C)-B \left (c^2-d^2\right )\right )+2 a b \left (-A \left (c^2-d^2\right )-2 B c d+c^2 C-C d^2\right )-b^2 \left (2 c d (A-C)-B \left (c^2-d^2\right )\right )\right )}{f \left (c^2+d^2\right )^2}-\frac {x \left (a^2 \left (-A \left (c^2-d^2\right )-2 B c d+c^2 C-C d^2\right )-2 a b \left (2 c d (A-C)-B \left (c^2-d^2\right )\right )-b^2 \left (-A \left (c^2-d^2\right )-2 B c d+c^2 C-C d^2\right )\right )}{\left (c^2+d^2\right )^2}-\frac {\left (A d^2-B c d+c^2 C\right ) (a+b \tan (e+f x))^2}{d f \left (c^2+d^2\right ) (c+d \tan (e+f x))}-\frac {(b c-a d) \left (a d^2 \left (2 c d (A-C)-B \left (c^2-d^2\right )\right )+b \left (2 A d^4-B c^3 d-3 B c d^3+2 c^4 C+4 c^2 C d^2\right )\right ) \log (c+d \tan (e+f x))}{d^3 f \left (c^2+d^2\right )^2}+\frac {b^2 \tan (e+f x) \left (d^2 (A+C)-B c d+2 c^2 C\right )}{d^2 f \left (c^2+d^2\right )} \]

[In]

Int[((a + b*Tan[e + f*x])^2*(A + B*Tan[e + f*x] + C*Tan[e + f*x]^2))/(c + d*Tan[e + f*x])^2,x]

[Out]

-(((a^2*(c^2*C - 2*B*c*d - C*d^2 - A*(c^2 - d^2)) - b^2*(c^2*C - 2*B*c*d - C*d^2 - A*(c^2 - d^2)) - 2*a*b*(2*c
*(A - C)*d - B*(c^2 - d^2)))*x)/(c^2 + d^2)^2) + ((2*a*b*(c^2*C - 2*B*c*d - C*d^2 - A*(c^2 - d^2)) + a^2*(2*c*
(A - C)*d - B*(c^2 - d^2)) - b^2*(2*c*(A - C)*d - B*(c^2 - d^2)))*Log[Cos[e + f*x]])/((c^2 + d^2)^2*f) - ((b*c
- a*d)*(b*(2*c^4*C - B*c^3*d + 4*c^2*C*d^2 - 3*B*c*d^3 + 2*A*d^4) + a*d^2*(2*c*(A - C)*d - B*(c^2 - d^2)))*Lo
g[c + d*Tan[e + f*x]])/(d^3*(c^2 + d^2)^2*f) + (b^2*(2*c^2*C - B*c*d + (A + C)*d^2)*Tan[e + f*x])/(d^2*(c^2 +
d^2)*f) - ((c^2*C - B*c*d + A*d^2)*(a + b*Tan[e + f*x])^2)/(d*(c^2 + d^2)*f*(c + d*Tan[e + f*x]))

Rule 31

Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]

Rule 3556

Int[tan[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Log[RemoveContent[Cos[c + d*x], x]]/d, x] /; FreeQ[{c, d}, x]

Rule 3698

Int[((a_.) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(m_.)*((A_) + (C_.)*tan[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Dist[
A/(b*f), Subst[Int[(a + x)^m, x], x, b*Tan[e + f*x]], x] /; FreeQ[{a, b, e, f, A, C, m}, x] && EqQ[A, C]

Rule 3707

Int[((A_) + (B_.)*tan[(e_.) + (f_.)*(x_)] + (C_.)*tan[(e_.) + (f_.)*(x_)]^2)/((a_.) + (b_.)*tan[(e_.) + (f_.)*
(x_)]), x_Symbol] :> Simp[(a*A + b*B - a*C)*(x/(a^2 + b^2)), x] + (Dist[(A*b^2 - a*b*B + a^2*C)/(a^2 + b^2), I
nt[(1 + Tan[e + f*x]^2)/(a + b*Tan[e + f*x]), x], x] - Dist[(A*b - a*B - b*C)/(a^2 + b^2), Int[Tan[e + f*x], x
], x]) /; FreeQ[{a, b, e, f, A, B, C}, x] && NeQ[A*b^2 - a*b*B + a^2*C, 0] && NeQ[a^2 + b^2, 0] && NeQ[A*b - a
*B - b*C, 0]

Rule 3718

Int[((a_) + (b_.)*tan[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*tan[(e_.) + (f_.)*(x_)])^(n_.)*((A_.) + (B_.)*tan[(e
_.) + (f_.)*(x_)] + (C_.)*tan[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[b*C*Tan[e + f*x]*((c + d*Tan[e + f*x])
^(n + 1)/(d*f*(n + 2))), x] - Dist[1/(d*(n + 2)), Int[(c + d*Tan[e + f*x])^n*Simp[b*c*C - a*A*d*(n + 2) - (A*b
+ a*B - b*C)*d*(n + 2)*Tan[e + f*x] - (a*C*d*(n + 2) - b*(c*C - B*d*(n + 2)))*Tan[e + f*x]^2, x], x], x] /; F
reeQ[{a, b, c, d, e, f, A, B, C, n}, x] && NeQ[b*c - a*d, 0] && NeQ[c^2 + d^2, 0] &&  !LtQ[n, -1]

Rule 3726

Int[((a_.) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*tan[(e_.) + (f_.)*(x_)])^(n_)*((A_.) + (B_.)*t
an[(e_.) + (f_.)*(x_)] + (C_.)*tan[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[(A*d^2 + c*(c*C - B*d))*(a + b*Ta
n[e + f*x])^m*((c + d*Tan[e + f*x])^(n + 1)/(d*f*(n + 1)*(c^2 + d^2))), x] - Dist[1/(d*(n + 1)*(c^2 + d^2)), I
nt[(a + b*Tan[e + f*x])^(m - 1)*(c + d*Tan[e + f*x])^(n + 1)*Simp[A*d*(b*d*m - a*c*(n + 1)) + (c*C - B*d)*(b*c
*m + a*d*(n + 1)) - d*(n + 1)*((A - C)*(b*c - a*d) + B*(a*c + b*d))*Tan[e + f*x] - b*(d*(B*c - A*d)*(m + n + 1
) - C*(c^2*m - d^2*(n + 1)))*Tan[e + f*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c -
a*d, 0] && NeQ[a^2 + b^2, 0] && NeQ[c^2 + d^2, 0] && GtQ[m, 0] && LtQ[n, -1]

Rubi steps \begin{align*} \text {integral}& = -\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))}+\frac {\int \frac {(a+b \tan (e+f x)) \left (A d (a c+2 b d)+(2 b c-a d) (c C-B d)+d ((A-C) (b c-a d)+B (a c+b d)) \tan (e+f x)+b \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan ^2(e+f x)\right )}{c+d \tan (e+f x)} \, dx}{d \left (c^2+d^2\right )} \\ & = \frac {b^2 \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan (e+f x)}{d^2 \left (c^2+d^2\right ) f}-\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))}-\frac {\int \frac {b^2 c \left (2 c^2 C-B c d+(A+C) d^2\right )-a d (A d (a c+2 b d)+(2 b c-a d) (c C-B d))-d^2 \left (2 a b (A c-c C+B d)+a^2 (B c-(A-C) d)-b^2 (B c-(A-C) d)\right ) \tan (e+f x)+b (2 b c C-b B d-2 a C d) \left (c^2+d^2\right ) \tan ^2(e+f x)}{c+d \tan (e+f x)} \, dx}{d^2 \left (c^2+d^2\right )} \\ & = -\frac {\left (a^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-b^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-2 a b \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) x}{\left (c^2+d^2\right )^2}+\frac {b^2 \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan (e+f x)}{d^2 \left (c^2+d^2\right ) f}-\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))}-\frac {\left (2 a b \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )+a^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )-b^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \int \tan (e+f x) \, dx}{\left (c^2+d^2\right )^2}-\frac {\left ((b c-a d) \left (b \left (2 c^4 C-B c^3 d+4 c^2 C d^2-3 B c d^3+2 A d^4\right )+a d^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right )\right ) \int \frac {1+\tan ^2(e+f x)}{c+d \tan (e+f x)} \, dx}{d^2 \left (c^2+d^2\right )^2} \\ & = -\frac {\left (a^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-b^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-2 a b \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) x}{\left (c^2+d^2\right )^2}+\frac {\left (2 a b \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )+a^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )-b^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \log (\cos (e+f x))}{\left (c^2+d^2\right )^2 f}+\frac {b^2 \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan (e+f x)}{d^2 \left (c^2+d^2\right ) f}-\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))}-\frac {\left ((b c-a d) \left (b \left (2 c^4 C-B c^3 d+4 c^2 C d^2-3 B c d^3+2 A d^4\right )+a d^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right )\right ) \text {Subst}\left (\int \frac {1}{c+x} \, dx,x,d \tan (e+f x)\right )}{d^3 \left (c^2+d^2\right )^2 f} \\ & = -\frac {\left (a^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-b^2 \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )-2 a b \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) x}{\left (c^2+d^2\right )^2}+\frac {\left (2 a b \left (c^2 C-2 B c d-C d^2-A \left (c^2-d^2\right )\right )+a^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )-b^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \log (\cos (e+f x))}{\left (c^2+d^2\right )^2 f}-\frac {(b c-a d) \left (b \left (2 c^4 C-B c^3 d+4 c^2 C d^2-3 B c d^3+2 A d^4\right )+a d^2 \left (2 c (A-C) d-B \left (c^2-d^2\right )\right )\right ) \log (c+d \tan (e+f x))}{d^3 \left (c^2+d^2\right )^2 f}+\frac {b^2 \left (2 c^2 C-B c d+(A+C) d^2\right ) \tan (e+f x)}{d^2 \left (c^2+d^2\right ) f}-\frac {\left (c^2 C-B c d+A d^2\right ) (a+b \tan (e+f x))^2}{d \left (c^2+d^2\right ) f (c+d \tan (e+f x))} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 5.57 (sec) , antiderivative size = 277, normalized size of antiderivative = 0.66 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\frac {\frac {(a+i b)^2 (-i A+B+i C) \log (i-\tan (e+f x))}{(c+i d)^2}+\frac {(a-i b)^2 (i A+B-i C) \log (i+\tan (e+f x))}{(c-i d)^2}+\frac {2 (-b c+a d) \left (b \left (2 c^4 C-B c^3 d+4 c^2 C d^2-3 B c d^3+2 A d^4\right )+a d^2 \left (2 c (A-C) d+B \left (-c^2+d^2\right )\right )\right ) \log (c+d \tan (e+f x))}{d^3 \left (c^2+d^2\right )^2}-\frac {2 (b c-a d)^2 \left (2 c^2 C-B c d+(A+C) d^2\right )}{d^3 \left (c^2+d^2\right ) (c+d \tan (e+f x))}+\frac {2 C (a+b \tan (e+f x))^2}{d (c+d \tan (e+f x))}}{2 f} \]

[In]

Integrate[((a + b*Tan[e + f*x])^2*(A + B*Tan[e + f*x] + C*Tan[e + f*x]^2))/(c + d*Tan[e + f*x])^2,x]

[Out]

(((a + I*b)^2*((-I)*A + B + I*C)*Log[I - Tan[e + f*x]])/(c + I*d)^2 + ((a - I*b)^2*(I*A + B - I*C)*Log[I + Tan
[e + f*x]])/(c - I*d)^2 + (2*(-(b*c) + a*d)*(b*(2*c^4*C - B*c^3*d + 4*c^2*C*d^2 - 3*B*c*d^3 + 2*A*d^4) + a*d^2
*(2*c*(A - C)*d + B*(-c^2 + d^2)))*Log[c + d*Tan[e + f*x]])/(d^3*(c^2 + d^2)^2) - (2*(b*c - a*d)^2*(2*c^2*C -
B*c*d + (A + C)*d^2))/(d^3*(c^2 + d^2)*(c + d*Tan[e + f*x])) + (2*C*(a + b*Tan[e + f*x])^2)/(d*(c + d*Tan[e +
f*x])))/(2*f)

Maple [A] (verified)

Time = 0.27 (sec) , antiderivative size = 552, normalized size of antiderivative = 1.32

method	result	size
derivativedivides	\(\frac {\frac {\tan \left (f x +e \right ) C \,b^{2}}{d^{2}}+\frac {\frac {\left (-2 A \,a^{2} c d +2 A a b \,c^{2}-2 A a b \,d^{2}+2 A \,b^{2} c d +B \,a^{2} c^{2}-B \,a^{2} d^{2}+4 B a b c d -B \,b^{2} c^{2}+B \,b^{2} d^{2}+2 C \,a^{2} c d -2 C a b \,c^{2}+2 C a b \,d^{2}-2 C \,b^{2} c d \right ) \ln \left (1+\tan \left (f x +e \right )^{2}\right )}{2}+\left (A \,a^{2} c^{2}-A \,a^{2} d^{2}+4 A a b c d -A \,b^{2} c^{2}+A \,b^{2} d^{2}+2 B \,a^{2} c d -2 B a b \,c^{2}+2 B a b \,d^{2}-2 B \,b^{2} c d -C \,a^{2} c^{2}+a^{2} C \,d^{2}-4 C a b c d +C \,b^{2} c^{2}-C \,b^{2} d^{2}\right ) \arctan \left (\tan \left (f x +e \right )\right )}{\left (c^{2}+d^{2}\right )^{2}}-\frac {A \,a^{2} d^{4}-2 A a b c \,d^{3}+A \,b^{2} c^{2} d^{2}-B \,a^{2} c \,d^{3}+2 B a b \,c^{2} d^{2}-B \,b^{2} c^{3} d +C \,a^{2} c^{2} d^{2}-2 C a b \,c^{3} d +C \,b^{2} c^{4}}{d^{3} \left (c^{2}+d^{2}\right ) \left (c +d \tan \left (f x +e \right )\right )}+\frac {\left (2 A \,a^{2} c \,d^{4}-2 A a b \,c^{2} d^{3}+2 A a b \,d^{5}-2 A \,b^{2} c \,d^{4}-B \,a^{2} c^{2} d^{3}+B \,a^{2} d^{5}-4 B a b c \,d^{4}+B \,b^{2} c^{4} d +3 B \,b^{2} c^{2} d^{3}-2 C \,a^{2} c \,d^{4}+2 C a b \,c^{4} d +6 C a b \,c^{2} d^{3}-2 C \,b^{2} c^{5}-4 C \,b^{2} c^{3} d^{2}\right ) \ln \left (c +d \tan \left (f x +e \right )\right )}{d^{3} \left (c^{2}+d^{2}\right )^{2}}}{f}\)	\(552\)
default	\(\frac {\frac {\tan \left (f x +e \right ) C \,b^{2}}{d^{2}}+\frac {\frac {\left (-2 A \,a^{2} c d +2 A a b \,c^{2}-2 A a b \,d^{2}+2 A \,b^{2} c d +B \,a^{2} c^{2}-B \,a^{2} d^{2}+4 B a b c d -B \,b^{2} c^{2}+B \,b^{2} d^{2}+2 C \,a^{2} c d -2 C a b \,c^{2}+2 C a b \,d^{2}-2 C \,b^{2} c d \right ) \ln \left (1+\tan \left (f x +e \right )^{2}\right )}{2}+\left (A \,a^{2} c^{2}-A \,a^{2} d^{2}+4 A a b c d -A \,b^{2} c^{2}+A \,b^{2} d^{2}+2 B \,a^{2} c d -2 B a b \,c^{2}+2 B a b \,d^{2}-2 B \,b^{2} c d -C \,a^{2} c^{2}+a^{2} C \,d^{2}-4 C a b c d +C \,b^{2} c^{2}-C \,b^{2} d^{2}\right ) \arctan \left (\tan \left (f x +e \right )\right )}{\left (c^{2}+d^{2}\right )^{2}}-\frac {A \,a^{2} d^{4}-2 A a b c \,d^{3}+A \,b^{2} c^{2} d^{2}-B \,a^{2} c \,d^{3}+2 B a b \,c^{2} d^{2}-B \,b^{2} c^{3} d +C \,a^{2} c^{2} d^{2}-2 C a b \,c^{3} d +C \,b^{2} c^{4}}{d^{3} \left (c^{2}+d^{2}\right ) \left (c +d \tan \left (f x +e \right )\right )}+\frac {\left (2 A \,a^{2} c \,d^{4}-2 A a b \,c^{2} d^{3}+2 A a b \,d^{5}-2 A \,b^{2} c \,d^{4}-B \,a^{2} c^{2} d^{3}+B \,a^{2} d^{5}-4 B a b c \,d^{4}+B \,b^{2} c^{4} d +3 B \,b^{2} c^{2} d^{3}-2 C \,a^{2} c \,d^{4}+2 C a b \,c^{4} d +6 C a b \,c^{2} d^{3}-2 C \,b^{2} c^{5}-4 C \,b^{2} c^{3} d^{2}\right ) \ln \left (c +d \tan \left (f x +e \right )\right )}{d^{3} \left (c^{2}+d^{2}\right )^{2}}}{f}\)	\(552\)
norman	\(\frac {\frac {c \left (A \,a^{2} c^{2}-A \,a^{2} d^{2}+4 A a b c d -A \,b^{2} c^{2}+A \,b^{2} d^{2}+2 B \,a^{2} c d -2 B a b \,c^{2}+2 B a b \,d^{2}-2 B \,b^{2} c d -C \,a^{2} c^{2}+a^{2} C \,d^{2}-4 C a b c d +C \,b^{2} c^{2}-C \,b^{2} d^{2}\right ) x}{c^{4}+2 c^{2} d^{2}+d^{4}}+\frac {C \,b^{2} \tan \left (f x +e \right )^{2}}{d f}+\frac {d \left (A \,a^{2} c^{2}-A \,a^{2} d^{2}+4 A a b c d -A \,b^{2} c^{2}+A \,b^{2} d^{2}+2 B \,a^{2} c d -2 B a b \,c^{2}+2 B a b \,d^{2}-2 B \,b^{2} c d -C \,a^{2} c^{2}+a^{2} C \,d^{2}-4 C a b c d +C \,b^{2} c^{2}-C \,b^{2} d^{2}\right ) x \tan \left (f x +e \right )}{c^{4}+2 c^{2} d^{2}+d^{4}}-\frac {A \,a^{2} d^{4}-2 A a b c \,d^{3}+A \,b^{2} c^{2} d^{2}-B \,a^{2} c \,d^{3}+2 B a b \,c^{2} d^{2}-B \,b^{2} c^{3} d +C \,a^{2} c^{2} d^{2}-2 C a b \,c^{3} d +2 C \,b^{2} c^{4}+C \,b^{2} c^{2} d^{2}}{f \,d^{3} \left (c^{2}+d^{2}\right )}}{c +d \tan \left (f x +e \right )}+\frac {\left (2 A \,a^{2} c \,d^{4}-2 A a b \,c^{2} d^{3}+2 A a b \,d^{5}-2 A \,b^{2} c \,d^{4}-B \,a^{2} c^{2} d^{3}+B \,a^{2} d^{5}-4 B a b c \,d^{4}+B \,b^{2} c^{4} d +3 B \,b^{2} c^{2} d^{3}-2 C \,a^{2} c \,d^{4}+2 C a b \,c^{4} d +6 C a b \,c^{2} d^{3}-2 C \,b^{2} c^{5}-4 C \,b^{2} c^{3} d^{2}\right ) \ln \left (c +d \tan \left (f x +e \right )\right )}{\left (c^{4}+2 c^{2} d^{2}+d^{4}\right ) d^{3} f}-\frac {\left (2 A \,a^{2} c d -2 A a b \,c^{2}+2 A a b \,d^{2}-2 A \,b^{2} c d -B \,a^{2} c^{2}+B \,a^{2} d^{2}-4 B a b c d +B \,b^{2} c^{2}-B \,b^{2} d^{2}-2 C \,a^{2} c d +2 C a b \,c^{2}-2 C a b \,d^{2}+2 C \,b^{2} c d \right ) \ln \left (1+\tan \left (f x +e \right )^{2}\right )}{2 f \left (c^{4}+2 c^{2} d^{2}+d^{4}\right )}\)	\(743\)
parallelrisch	\(\text {Expression too large to display}\)	\(2062\)
risch	\(\text {Expression too large to display}\)	\(2406\)

[In]

int((a+b*tan(f*x+e))^2*(A+B*tan(f*x+e)+C*tan(f*x+e)^2)/(c+d*tan(f*x+e))^2,x,method=_RETURNVERBOSE)

[Out]

1/f*(tan(f*x+e)*C*b^2/d^2+1/(c^2+d^2)^2*(1/2*(-2*A*a^2*c*d+2*A*a*b*c^2-2*A*a*b*d^2+2*A*b^2*c*d+B*a^2*c^2-B*a^2
*d^2+4*B*a*b*c*d-B*b^2*c^2+B*b^2*d^2+2*C*a^2*c*d-2*C*a*b*c^2+2*C*a*b*d^2-2*C*b^2*c*d)*ln(1+tan(f*x+e)^2)+(A*a^
2*c^2-A*a^2*d^2+4*A*a*b*c*d-A*b^2*c^2+A*b^2*d^2+2*B*a^2*c*d-2*B*a*b*c^2+2*B*a*b*d^2-2*B*b^2*c*d-C*a^2*c^2+C*a^
2*d^2-4*C*a*b*c*d+C*b^2*c^2-C*b^2*d^2)*arctan(tan(f*x+e)))-1/d^3*(A*a^2*d^4-2*A*a*b*c*d^3+A*b^2*c^2*d^2-B*a^2*
c*d^3+2*B*a*b*c^2*d^2-B*b^2*c^3*d+C*a^2*c^2*d^2-2*C*a*b*c^3*d+C*b^2*c^4)/(c^2+d^2)/(c+d*tan(f*x+e))+1/d^3*(2*A
*a^2*c*d^4-2*A*a*b*c^2*d^3+2*A*a*b*d^5-2*A*b^2*c*d^4-B*a^2*c^2*d^3+B*a^2*d^5-4*B*a*b*c*d^4+B*b^2*c^4*d+3*B*b^2
*c^2*d^3-2*C*a^2*c*d^4+2*C*a*b*c^4*d+6*C*a*b*c^2*d^3-2*C*b^2*c^5-4*C*b^2*c^3*d^2)/(c^2+d^2)^2*ln(c+d*tan(f*x+e
)))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 939 vs. \(2 (417) = 834\).

Time = 0.56 (sec) , antiderivative size = 939, normalized size of antiderivative = 2.25 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=-\frac {2 \, C b^{2} c^{4} d^{2} + 2 \, A a^{2} d^{6} - 2 \, {\left (2 \, C a b + B b^{2}\right )} c^{3} d^{3} + 2 \, {\left (C a^{2} + 2 \, B a b + A b^{2}\right )} c^{2} d^{4} - 2 \, {\left (B a^{2} + 2 \, A a b\right )} c d^{5} - 2 \, {\left ({\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} c^{3} d^{3} + 2 \, {\left (B a^{2} + 2 \, {\left (A - C\right )} a b - B b^{2}\right )} c^{2} d^{4} - {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} c d^{5}\right )} f x - 2 \, {\left (C b^{2} c^{4} d^{2} + 2 \, C b^{2} c^{2} d^{4} + C b^{2} d^{6}\right )} \tan \left (f x + e\right )^{2} + {\left (2 \, C b^{2} c^{6} + 4 \, C b^{2} c^{4} d^{2} - {\left (2 \, C a b + B b^{2}\right )} c^{5} d + {\left (B a^{2} + 2 \, {\left (A - 3 \, C\right )} a b - 3 \, B b^{2}\right )} c^{3} d^{3} - 2 \, {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - A b^{2}\right )} c^{2} d^{4} - {\left (B a^{2} + 2 \, A a b\right )} c d^{5} + {\left (2 \, C b^{2} c^{5} d + 4 \, C b^{2} c^{3} d^{3} - {\left (2 \, C a b + B b^{2}\right )} c^{4} d^{2} + {\left (B a^{2} + 2 \, {\left (A - 3 \, C\right )} a b - 3 \, B b^{2}\right )} c^{2} d^{4} - 2 \, {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - A b^{2}\right )} c d^{5} - {\left (B a^{2} + 2 \, A a b\right )} d^{6}\right )} \tan \left (f x + e\right )\right )} \log \left (\frac {d^{2} \tan \left (f x + e\right )^{2} + 2 \, c d \tan \left (f x + e\right ) + c^{2}}{\tan \left (f x + e\right )^{2} + 1}\right ) - {\left (2 \, C b^{2} c^{6} + 4 \, C b^{2} c^{4} d^{2} + 2 \, C b^{2} c^{2} d^{4} - {\left (2 \, C a b + B b^{2}\right )} c^{5} d - 2 \, {\left (2 \, C a b + B b^{2}\right )} c^{3} d^{3} - {\left (2 \, C a b + B b^{2}\right )} c d^{5} + {\left (2 \, C b^{2} c^{5} d + 4 \, C b^{2} c^{3} d^{3} + 2 \, C b^{2} c d^{5} - {\left (2 \, C a b + B b^{2}\right )} c^{4} d^{2} - 2 \, {\left (2 \, C a b + B b^{2}\right )} c^{2} d^{4} - {\left (2 \, C a b + B b^{2}\right )} d^{6}\right )} \tan \left (f x + e\right )\right )} \log \left (\frac {1}{\tan \left (f x + e\right )^{2} + 1}\right ) - 2 \, {\left (2 \, C b^{2} c^{5} d - {\left (2 \, C a b + B b^{2}\right )} c^{4} d^{2} + {\left (C a^{2} + 2 \, B a b + {\left (A + 2 \, C\right )} b^{2}\right )} c^{3} d^{3} - {\left (B a^{2} + 2 \, A a b\right )} c^{2} d^{4} + {\left (A a^{2} + C b^{2}\right )} c d^{5} + {\left ({\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} c^{2} d^{4} + 2 \, {\left (B a^{2} + 2 \, {\left (A - C\right )} a b - B b^{2}\right )} c d^{5} - {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} d^{6}\right )} f x\right )} \tan \left (f x + e\right )}{2 \, {\left ({\left (c^{4} d^{4} + 2 \, c^{2} d^{6} + d^{8}\right )} f \tan \left (f x + e\right ) + {\left (c^{5} d^{3} + 2 \, c^{3} d^{5} + c d^{7}\right )} f\right )}} \]

[In]

integrate((a+b*tan(f*x+e))^2*(A+B*tan(f*x+e)+C*tan(f*x+e)^2)/(c+d*tan(f*x+e))^2,x, algorithm="fricas")

[Out]

-1/2*(2*C*b^2*c^4*d^2 + 2*A*a^2*d^6 - 2*(2*C*a*b + B*b^2)*c^3*d^3 + 2*(C*a^2 + 2*B*a*b + A*b^2)*c^2*d^4 - 2*(B
*a^2 + 2*A*a*b)*c*d^5 - 2*(((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*c^3*d^3 + 2*(B*a^2 + 2*(A - C)*a*b - B*b^2)*c
^2*d^4 - ((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*c*d^5)*f*x - 2*(C*b^2*c^4*d^2 + 2*C*b^2*c^2*d^4 + C*b^2*d^6)*ta
n(f*x + e)^2 + (2*C*b^2*c^6 + 4*C*b^2*c^4*d^2 - (2*C*a*b + B*b^2)*c^5*d + (B*a^2 + 2*(A - 3*C)*a*b - 3*B*b^2)*
c^3*d^3 - 2*((A - C)*a^2 - 2*B*a*b - A*b^2)*c^2*d^4 - (B*a^2 + 2*A*a*b)*c*d^5 + (2*C*b^2*c^5*d + 4*C*b^2*c^3*d
^3 - (2*C*a*b + B*b^2)*c^4*d^2 + (B*a^2 + 2*(A - 3*C)*a*b - 3*B*b^2)*c^2*d^4 - 2*((A - C)*a^2 - 2*B*a*b - A*b^
2)*c*d^5 - (B*a^2 + 2*A*a*b)*d^6)*tan(f*x + e))*log((d^2*tan(f*x + e)^2 + 2*c*d*tan(f*x + e) + c^2)/(tan(f*x +
e)^2 + 1)) - (2*C*b^2*c^6 + 4*C*b^2*c^4*d^2 + 2*C*b^2*c^2*d^4 - (2*C*a*b + B*b^2)*c^5*d - 2*(2*C*a*b + B*b^2)
*c^3*d^3 - (2*C*a*b + B*b^2)*c*d^5 + (2*C*b^2*c^5*d + 4*C*b^2*c^3*d^3 + 2*C*b^2*c*d^5 - (2*C*a*b + B*b^2)*c^4*
d^2 - 2*(2*C*a*b + B*b^2)*c^2*d^4 - (2*C*a*b + B*b^2)*d^6)*tan(f*x + e))*log(1/(tan(f*x + e)^2 + 1)) - 2*(2*C*
b^2*c^5*d - (2*C*a*b + B*b^2)*c^4*d^2 + (C*a^2 + 2*B*a*b + (A + 2*C)*b^2)*c^3*d^3 - (B*a^2 + 2*A*a*b)*c^2*d^4
+ (A*a^2 + C*b^2)*c*d^5 + (((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*c^2*d^4 + 2*(B*a^2 + 2*(A - C)*a*b - B*b^2)*c
*d^5 - ((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*d^6)*f*x)*tan(f*x + e))/((c^4*d^4 + 2*c^2*d^6 + d^8)*f*tan(f*x +
e) + (c^5*d^3 + 2*c^3*d^5 + c*d^7)*f)

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 3.01 (sec) , antiderivative size = 16225, normalized size of antiderivative = 38.91 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\text {Too large to display} \]

[In]

integrate((a+b*tan(f*x+e))**2*(A+B*tan(f*x+e)+C*tan(f*x+e)**2)/(c+d*tan(f*x+e))**2,x)

[Out]

Piecewise((zoo*x*(a + b*tan(e))**2*(A + B*tan(e) + C*tan(e)**2)/tan(e)**2, Eq(c, 0) & Eq(d, 0) & Eq(f, 0)), ((
A*a**2*x + A*a*b*log(tan(e + f*x)**2 + 1)/f - A*b**2*x + A*b**2*tan(e + f*x)/f + B*a**2*log(tan(e + f*x)**2 +
1)/(2*f) - 2*B*a*b*x + 2*B*a*b*tan(e + f*x)/f - B*b**2*log(tan(e + f*x)**2 + 1)/(2*f) + B*b**2*tan(e + f*x)**2
/(2*f) - C*a**2*x + C*a**2*tan(e + f*x)/f - C*a*b*log(tan(e + f*x)**2 + 1)/f + C*a*b*tan(e + f*x)**2/f + C*b**
2*x + C*b**2*tan(e + f*x)**3/(3*f) - C*b**2*tan(e + f*x)/f)/c**2, Eq(d, 0)), (-A*a**2*f*x*tan(e + f*x)**2/(4*d
**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*I*A*a**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*
x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + A*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) -
4*d**2*f) - A*a**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*I*A*a**2/
(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*I*A*a*b*f*x*tan(e + f*x)**2/(4*d**2*f*tan(
e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 4*A*a*b*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d
**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*a*b*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f)
+ 2*I*A*a*b*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + A*b**2*f*x*tan(e +
f*x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*b**2*f*x*tan(e + f*x)/(4*d**2
*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - A*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*t
an(e + f*x) - 4*d**2*f) - 3*A*b**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f
) + 2*I*A*b**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + I*B*a**2*f*x*tan(e + f*x)**2/
(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*B*a**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e +
f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - I*B*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*
x) - 4*d**2*f) + I*B*a**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*B*a
*b*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 4*I*B*a*b*f*x*tan(e +
f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*B*a*b*f*x/(4*d**2*f*tan(e + f*x)**2
- 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 6*B*a*b*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*
x) - 4*d**2*f) + 4*I*B*a*b/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 3*I*B*b**2*f*x*ta
n(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 6*B*b**2*f*x*tan(e + f*x)/(4*d
**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 3*I*B*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d
**2*f*tan(e + f*x) - 4*d**2*f) + 2*B*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 -
8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 4*I*B*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)
**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*B*b**2*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 - 8*I*
d**2*f*tan(e + f*x) - 4*d**2*f) - 5*I*B*b**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x)
- 4*d**2*f) - 4*B*b**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + C*a**2*f*x*tan(e + f*
x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*C*a**2*f*x*tan(e + f*x)/(4*d**2*f*
tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - C*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(
e + f*x) - 4*d**2*f) - 3*C*a**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) +
2*I*C*a**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 6*I*C*a*b*f*x*tan(e + f*x)**2/(4
*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 12*C*a*b*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*
x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 6*I*C*a*b*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x
) - 4*d**2*f) + 4*C*a*b*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e
+ f*x) - 4*d**2*f) - 8*I*C*a*b*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*ta
n(e + f*x) - 4*d**2*f) - 4*C*a*b*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x)
- 4*d**2*f) - 10*I*C*a*b*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 8*C*a*
b/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 9*C*b**2*f*x*tan(e + f*x)**2/(4*d**2*f*tan
(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 18*I*C*b**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 -
8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 9*C*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**
2*f) + 4*I*C*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x)
- 4*d**2*f) + 8*C*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f
*x) - 4*d**2*f) - 4*I*C*b**2*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*
d**2*f) + 4*C*b**2*tan(e + f*x)**3/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 19*C*b**2
*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 14*I*C*b**2/(4*d**2*f*tan(e +
f*x)**2 - 8*I*d**2*f*tan(e + f*x) - 4*d**2*f), Eq(c, -I*d)), (-A*a**2*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e + f*
x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*a**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2
*f*tan(e + f*x) - 4*d**2*f) + A*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - A*a
**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*a**2/(4*d**2*f*tan(e
+ f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*a*b*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8*
I*d**2*f*tan(e + f*x) - 4*d**2*f) + 4*A*a*b*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*
x) - 4*d**2*f) + 2*I*A*a*b*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*a*b*tan
(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + A*b**2*f*x*tan(e + f*x)**2/(4*d**2
*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*I*A*b**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)*
*2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - A*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*
d**2*f) - 3*A*b**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*A*b**2/(
4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - I*B*a**2*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e
+ f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*B*a**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d*
*2*f*tan(e + f*x) - 4*d**2*f) + I*B*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) -
I*B*a**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*B*a*b*f*x*tan(e + f
*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 4*I*B*a*b*f*x*tan(e + f*x)/(4*d**2*f*
tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*B*a*b*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan
(e + f*x) - 4*d**2*f) - 6*B*a*b*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) -
4*I*B*a*b/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 3*I*B*b**2*f*x*tan(e + f*x)**2/(4
*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 6*B*b**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*
x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 3*I*B*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*
x) - 4*d**2*f) + 2*B*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(
e + f*x) - 4*d**2*f) + 4*I*B*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f
*tan(e + f*x) - 4*d**2*f) - 2*B*b**2*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f
*x) - 4*d**2*f) + 5*I*B*b**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 4*
B*b**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + C*a**2*f*x*tan(e + f*x)**2/(4*d**2*f*
tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 2*I*C*a**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2
+ 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - C*a**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**
2*f) - 3*C*a**2*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 2*I*C*a**2/(4*d
**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 6*I*C*a*b*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e +
f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 12*C*a*b*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2
*f*tan(e + f*x) - 4*d**2*f) + 6*I*C*a*b*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) +
4*C*a*b*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*
f) + 8*I*C*a*b*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d
**2*f) - 4*C*a*b*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 10
*I*C*a*b*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 8*C*a*b/(4*d**2*f*tan(
e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 9*C*b**2*f*x*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8
*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 18*I*C*b**2*f*x*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e
+ f*x) - 4*d**2*f) + 9*C*b**2*f*x/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) - 4*I*C*b**
2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 8
*C*b**2*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f)
+ 4*I*C*b**2*log(tan(e + f*x)**2 + 1)/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 4*C*b*
*2*tan(e + f*x)**3/(4*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 19*C*b**2*tan(e + f*x)/(4
*d**2*f*tan(e + f*x)**2 + 8*I*d**2*f*tan(e + f*x) - 4*d**2*f) + 14*I*C*b**2/(4*d**2*f*tan(e + f*x)**2 + 8*I*d*
*2*f*tan(e + f*x) - 4*d**2*f), Eq(c, I*d)), (x*(a + b*tan(e))**2*(A + B*tan(e) + C*tan(e)**2)/(c + d*tan(e))**
2, Eq(f, 0)), (2*A*a**2*c**3*d**3*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**
6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*a**2*c**2*d**4*f*x*tan(e + f*x)/(2*c**5*d**3*f +
2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x))
+ 4*A*a**2*c**2*d**4*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c
**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*a**2*c**2*d**4*log(tan(e + f*x)**2 + 1)/(2
*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f
*tan(e + f*x)) - 2*A*a**2*c**2*d**4/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*
f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*a**2*c*d**5*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(
e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*A*a**2*c*d**5*
log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6
*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*a**2*c*d**5*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/
(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8
*f*tan(e + f*x)) - 2*A*a**2*d**6*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f
+ 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*a**2*d**6/(2*c**5*d**3*f + 2*c**4*d**
4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*A*a*b*
c**3*d**3*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*
tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*a*b*c**3*d**3*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f
+ 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x
)) + 4*A*a*b*c**3*d**3/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x
) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 8*A*a*b*c**2*d**4*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) +
4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*A*a*b*c**2*d**4*log(c/d +
tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e
+ f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*a*b*c**2*d**4*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5
*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(
e + f*x)) + 8*A*a*b*c*d**5*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c*
*2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*A*a*b*c*d**5*log(c/d + tan(e + f*x))/(2*c**5*
d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e
+ f*x)) - 2*A*a*b*c*d**5*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f
+ 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*A*a*b*c*d**5/(2*c**5*d**3*f + 2*c**4*d
**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*A*a*
b*d**6*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c*
*2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*a*b*d**6*log(tan(e + f*x)**2 + 1)*tan(e + f
*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*
d**8*f*tan(e + f*x)) - 2*A*b**2*c**4*d**2/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2
*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*b**2*c**3*d**3*f*x/(2*c**5*d**3*f + 2*c**4*d*
*4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*b**
2*c**2*d**4*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e
+ f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*A*b**2*c**2*d**4*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*
c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) +
2*A*b**2*c**2*d**4*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c*
*2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*A*b**2*c**2*d**4/(2*c**5*d**3*f + 2*c**4*d**4
*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*b**2*
c*d**5*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7
*f + 2*d**8*f*tan(e + f*x)) - 4*A*b**2*c*d**5*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**
4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*b**2
*c*d**5*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*
c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*A*b**2*d**6*f*x*tan(e + f*x)/(2*c**5*d**3*f
+ 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x
)) - 2*B*a**2*c**3*d**3*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f +
4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + B*a**2*c**3*d**3*log(tan(e + f*x)**2 + 1)/(
2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*
f*tan(e + f*x)) + 2*B*a**2*c**3*d**3/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6
*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*B*a**2*c**2*d**4*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*
tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*B*a**2*c**
2*d**4*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c*
*2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + B*a**2*c**2*d**4*log(tan(e + f*x)**2 + 1)*tan(e
+ f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f
+ 2*d**8*f*tan(e + f*x)) + 4*B*a**2*c*d**5*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c*
*3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*B*a**2*c*d**5*log(c/d + tan(e
+ f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f
+ 2*d**8*f*tan(e + f*x)) - B*a**2*c*d**5*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x)
+ 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*B*a**2*c*d**5/(2*c**5*
d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e
+ f*x)) + 2*B*a**2*d**6*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*
c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - B*a**2*d**6*log(tan(e + f*x)*
*2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x)
+ 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*B*a*b*c**4*d**2/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3
*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*B*a*b*c**3*d**3*f*x/(2*c**5*d**
3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e +
f*x)) - 4*B*a*b*c**2*d**4*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**
2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 8*B*a*b*c**2*d**4*log(c/d + tan(e + f*x))/(2*c**
5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan
(e + f*x)) + 4*B*a*b*c**2*d**4*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d
**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*B*a*b*c**2*d**4/(2*c**5*d**3*f +
2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x))
+ 4*B*a*b*c*d**5*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x)
+ 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 8*B*a*b*c*d**5*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2
*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) +
4*B*a*b*c*d**5*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**
5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*B*a*b*d**6*f*x*tan(e + f*x)/(2*c**5
*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(
e + f*x)) + 2*B*b**2*c**5*d*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*
f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*B*b**2*c**5*d/(2*c**5*d**3*f + 2*c**4
*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*B*
b**2*c**4*d**2*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*
f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 6*B*b**2*c**3*d**3*log(c/d + tan(e + f*
x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*
d**8*f*tan(e + f*x)) - B*b**2*c**3*d**3*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) +
4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*B*b**2*c**3*d**3/(2*c**5
*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(
e + f*x)) - 4*B*b**2*c**2*d**4*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f
*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 6*B*b**2*c**2*d**4*log(c/d + tan(e + f*x))*tan(e + f*x)/
(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8
*f*tan(e + f*x)) - B*b**2*c**2*d**4*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e
+ f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*B*b**2*c*d**5*f
*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c
*d**7*f + 2*d**8*f*tan(e + f*x)) + B*b**2*c*d**5*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e
+ f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + B*b**2*d**6*log(t
an(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*t
an(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*a**2*c**4*d**2/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e +
f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*a**2*c**3*d**3*f
*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d
**8*f*tan(e + f*x)) - 2*C*a**2*c**2*d**4*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3
*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*C*a**2*c**2*d**4*log(c/d + tan(
e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*
f + 2*d**8*f*tan(e + f*x)) + 2*C*a**2*c**2*d**4*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e
+ f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*a**2*c**2*d**4
/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**
8*f*tan(e + f*x)) + 2*C*a**2*c*d**5*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d
**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*C*a**2*c*d**5*log(c/d + tan(e + f*x))*tan(e + f*x
)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d*
*8*f*tan(e + f*x)) + 2*C*a**2*c*d**5*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(
e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C*a**2*d**6*f*
x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*
d**7*f + 2*d**8*f*tan(e + f*x)) + 4*C*a*b*c**5*d*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e
+ f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*C*a*b*c**5*d/(2*
c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*
tan(e + f*x)) + 4*C*a*b*c**4*d**2*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e +
f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 12*C*a*b*c**3*d**3*l
og(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x
) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*a*b*c**3*d**3*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d
**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*C*a*
b*c**3*d**3/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**
7*f + 2*d**8*f*tan(e + f*x)) - 8*C*a*b*c**2*d**4*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5
*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 12*C*a*b*c**2*d**4*log(c/d + tan(e + f
*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2
*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*a*b*c**2*d**4*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f +
2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x))
- 8*C*a*b*c*d**5*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*
tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C*a*b*c*d**5*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f +
2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x))
+ 2*C*a*b*d**6*log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5
*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*C*b**2*c**6*log(c/d + tan(e + f*x))/
(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8
*f*tan(e + f*x)) - 4*C*b**2*c**6/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*t
an(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 4*C*b**2*c**5*d*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c*
*5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*ta
n(e + f*x)) - 8*C*b**2*c**4*d**2*log(c/d + tan(e + f*x))/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*
d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C*b**2*c**4*d**2*tan(e + f*x)**2
/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**
8*f*tan(e + f*x)) - 6*C*b**2*c**4*d**2/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d*
*6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C*b**2*c**3*d**3*f*x/(2*c**5*d**3*f + 2*c**4*d**4*
f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 8*C*b**2*c
**3*d**3*log(c/d + tan(e + f*x))*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*
c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C*b**2*c**2*d**4*f*x*tan(e + f*x)/(2*c**5*d
**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e
+ f*x)) - 2*C*b**2*c**2*d**4*log(tan(e + f*x)**2 + 1)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**
5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 4*C*b**2*c**2*d**4*tan(e + f*x)**2/(2
*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f
*tan(e + f*x)) - 2*C*b**2*c**2*d**4/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*
f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*b**2*c*d**5*f*x/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(
e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*b**2*c*d**5*
log(tan(e + f*x)**2 + 1)*tan(e + f*x)/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**
6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) - 2*C*b**2*d**6*f*x*tan(e + f*x)/(2*c**5*d**3*f + 2*c**
4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e + f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)) + 2*C
*b**2*d**6*tan(e + f*x)**2/(2*c**5*d**3*f + 2*c**4*d**4*f*tan(e + f*x) + 4*c**3*d**5*f + 4*c**2*d**6*f*tan(e +
f*x) + 2*c*d**7*f + 2*d**8*f*tan(e + f*x)), True))

Maxima [A] (verification not implemented)

none

Time = 0.34 (sec) , antiderivative size = 493, normalized size of antiderivative = 1.18 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\frac {\frac {2 \, C b^{2} \tan \left (f x + e\right )}{d^{2}} + \frac {2 \, {\left ({\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} c^{2} + 2 \, {\left (B a^{2} + 2 \, {\left (A - C\right )} a b - B b^{2}\right )} c d - {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} d^{2}\right )} {\left (f x + e\right )}}{c^{4} + 2 \, c^{2} d^{2} + d^{4}} - \frac {2 \, {\left (2 \, C b^{2} c^{5} + 4 \, C b^{2} c^{3} d^{2} - {\left (2 \, C a b + B b^{2}\right )} c^{4} d + {\left (B a^{2} + 2 \, {\left (A - 3 \, C\right )} a b - 3 \, B b^{2}\right )} c^{2} d^{3} - 2 \, {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - A b^{2}\right )} c d^{4} - {\left (B a^{2} + 2 \, A a b\right )} d^{5}\right )} \log \left (d \tan \left (f x + e\right ) + c\right )}{c^{4} d^{3} + 2 \, c^{2} d^{5} + d^{7}} + \frac {{\left ({\left (B a^{2} + 2 \, {\left (A - C\right )} a b - B b^{2}\right )} c^{2} - 2 \, {\left ({\left (A - C\right )} a^{2} - 2 \, B a b - {\left (A - C\right )} b^{2}\right )} c d - {\left (B a^{2} + 2 \, {\left (A - C\right )} a b - B b^{2}\right )} d^{2}\right )} \log \left (\tan \left (f x + e\right )^{2} + 1\right )}{c^{4} + 2 \, c^{2} d^{2} + d^{4}} - \frac {2 \, {\left (C b^{2} c^{4} + A a^{2} d^{4} - {\left (2 \, C a b + B b^{2}\right )} c^{3} d + {\left (C a^{2} + 2 \, B a b + A b^{2}\right )} c^{2} d^{2} - {\left (B a^{2} + 2 \, A a b\right )} c d^{3}\right )}}{c^{3} d^{3} + c d^{5} + {\left (c^{2} d^{4} + d^{6}\right )} \tan \left (f x + e\right )}}{2 \, f} \]

[In]

integrate((a+b*tan(f*x+e))^2*(A+B*tan(f*x+e)+C*tan(f*x+e)^2)/(c+d*tan(f*x+e))^2,x, algorithm="maxima")

[Out]

1/2*(2*C*b^2*tan(f*x + e)/d^2 + 2*(((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*c^2 + 2*(B*a^2 + 2*(A - C)*a*b - B*b^
2)*c*d - ((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*d^2)*(f*x + e)/(c^4 + 2*c^2*d^2 + d^4) - 2*(2*C*b^2*c^5 + 4*C*b
^2*c^3*d^2 - (2*C*a*b + B*b^2)*c^4*d + (B*a^2 + 2*(A - 3*C)*a*b - 3*B*b^2)*c^2*d^3 - 2*((A - C)*a^2 - 2*B*a*b
- A*b^2)*c*d^4 - (B*a^2 + 2*A*a*b)*d^5)*log(d*tan(f*x + e) + c)/(c^4*d^3 + 2*c^2*d^5 + d^7) + ((B*a^2 + 2*(A -
C)*a*b - B*b^2)*c^2 - 2*((A - C)*a^2 - 2*B*a*b - (A - C)*b^2)*c*d - (B*a^2 + 2*(A - C)*a*b - B*b^2)*d^2)*log(
tan(f*x + e)^2 + 1)/(c^4 + 2*c^2*d^2 + d^4) - 2*(C*b^2*c^4 + A*a^2*d^4 - (2*C*a*b + B*b^2)*c^3*d + (C*a^2 + 2*
B*a*b + A*b^2)*c^2*d^2 - (B*a^2 + 2*A*a*b)*c*d^3)/(c^3*d^3 + c*d^5 + (c^2*d^4 + d^6)*tan(f*x + e)))/f

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 893 vs. \(2 (417) = 834\).

Time = 0.78 (sec) , antiderivative size = 893, normalized size of antiderivative = 2.14 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\frac {\frac {2 \, C b^{2} \tan \left (f x + e\right )}{d^{2}} + \frac {2 \, {\left (A a^{2} c^{2} - C a^{2} c^{2} - 2 \, B a b c^{2} - A b^{2} c^{2} + C b^{2} c^{2} + 2 \, B a^{2} c d + 4 \, A a b c d - 4 \, C a b c d - 2 \, B b^{2} c d - A a^{2} d^{2} + C a^{2} d^{2} + 2 \, B a b d^{2} + A b^{2} d^{2} - C b^{2} d^{2}\right )} {\left (f x + e\right )}}{c^{4} + 2 \, c^{2} d^{2} + d^{4}} + \frac {{\left (B a^{2} c^{2} + 2 \, A a b c^{2} - 2 \, C a b c^{2} - B b^{2} c^{2} - 2 \, A a^{2} c d + 2 \, C a^{2} c d + 4 \, B a b c d + 2 \, A b^{2} c d - 2 \, C b^{2} c d - B a^{2} d^{2} - 2 \, A a b d^{2} + 2 \, C a b d^{2} + B b^{2} d^{2}\right )} \log \left (\tan \left (f x + e\right )^{2} + 1\right )}{c^{4} + 2 \, c^{2} d^{2} + d^{4}} - \frac {2 \, {\left (2 \, C b^{2} c^{5} - 2 \, C a b c^{4} d - B b^{2} c^{4} d + 4 \, C b^{2} c^{3} d^{2} + B a^{2} c^{2} d^{3} + 2 \, A a b c^{2} d^{3} - 6 \, C a b c^{2} d^{3} - 3 \, B b^{2} c^{2} d^{3} - 2 \, A a^{2} c d^{4} + 2 \, C a^{2} c d^{4} + 4 \, B a b c d^{4} + 2 \, A b^{2} c d^{4} - B a^{2} d^{5} - 2 \, A a b d^{5}\right )} \log \left ({\left | d \tan \left (f x + e\right ) + c \right |}\right )}{c^{4} d^{3} + 2 \, c^{2} d^{5} + d^{7}} + \frac {2 \, {\left (2 \, C b^{2} c^{5} d \tan \left (f x + e\right ) - 2 \, C a b c^{4} d^{2} \tan \left (f x + e\right ) - B b^{2} c^{4} d^{2} \tan \left (f x + e\right ) + 4 \, C b^{2} c^{3} d^{3} \tan \left (f x + e\right ) + B a^{2} c^{2} d^{4} \tan \left (f x + e\right ) + 2 \, A a b c^{2} d^{4} \tan \left (f x + e\right ) - 6 \, C a b c^{2} d^{4} \tan \left (f x + e\right ) - 3 \, B b^{2} c^{2} d^{4} \tan \left (f x + e\right ) - 2 \, A a^{2} c d^{5} \tan \left (f x + e\right ) + 2 \, C a^{2} c d^{5} \tan \left (f x + e\right ) + 4 \, B a b c d^{5} \tan \left (f x + e\right ) + 2 \, A b^{2} c d^{5} \tan \left (f x + e\right ) - B a^{2} d^{6} \tan \left (f x + e\right ) - 2 \, A a b d^{6} \tan \left (f x + e\right ) + C b^{2} c^{6} - C a^{2} c^{4} d^{2} - 2 \, B a b c^{4} d^{2} - A b^{2} c^{4} d^{2} + 3 \, C b^{2} c^{4} d^{2} + 2 \, B a^{2} c^{3} d^{3} + 4 \, A a b c^{3} d^{3} - 4 \, C a b c^{3} d^{3} - 2 \, B b^{2} c^{3} d^{3} - 3 \, A a^{2} c^{2} d^{4} + C a^{2} c^{2} d^{4} + 2 \, B a b c^{2} d^{4} + A b^{2} c^{2} d^{4} - A a^{2} d^{6}\right )}}{{\left (c^{4} d^{3} + 2 \, c^{2} d^{5} + d^{7}\right )} {\left (d \tan \left (f x + e\right ) + c\right )}}}{2 \, f} \]

[In]

integrate((a+b*tan(f*x+e))^2*(A+B*tan(f*x+e)+C*tan(f*x+e)^2)/(c+d*tan(f*x+e))^2,x, algorithm="giac")

[Out]

1/2*(2*C*b^2*tan(f*x + e)/d^2 + 2*(A*a^2*c^2 - C*a^2*c^2 - 2*B*a*b*c^2 - A*b^2*c^2 + C*b^2*c^2 + 2*B*a^2*c*d +
4*A*a*b*c*d - 4*C*a*b*c*d - 2*B*b^2*c*d - A*a^2*d^2 + C*a^2*d^2 + 2*B*a*b*d^2 + A*b^2*d^2 - C*b^2*d^2)*(f*x +
e)/(c^4 + 2*c^2*d^2 + d^4) + (B*a^2*c^2 + 2*A*a*b*c^2 - 2*C*a*b*c^2 - B*b^2*c^2 - 2*A*a^2*c*d + 2*C*a^2*c*d +
4*B*a*b*c*d + 2*A*b^2*c*d - 2*C*b^2*c*d - B*a^2*d^2 - 2*A*a*b*d^2 + 2*C*a*b*d^2 + B*b^2*d^2)*log(tan(f*x + e)
^2 + 1)/(c^4 + 2*c^2*d^2 + d^4) - 2*(2*C*b^2*c^5 - 2*C*a*b*c^4*d - B*b^2*c^4*d + 4*C*b^2*c^3*d^2 + B*a^2*c^2*d
^3 + 2*A*a*b*c^2*d^3 - 6*C*a*b*c^2*d^3 - 3*B*b^2*c^2*d^3 - 2*A*a^2*c*d^4 + 2*C*a^2*c*d^4 + 4*B*a*b*c*d^4 + 2*A
*b^2*c*d^4 - B*a^2*d^5 - 2*A*a*b*d^5)*log(abs(d*tan(f*x + e) + c))/(c^4*d^3 + 2*c^2*d^5 + d^7) + 2*(2*C*b^2*c^
5*d*tan(f*x + e) - 2*C*a*b*c^4*d^2*tan(f*x + e) - B*b^2*c^4*d^2*tan(f*x + e) + 4*C*b^2*c^3*d^3*tan(f*x + e) +
B*a^2*c^2*d^4*tan(f*x + e) + 2*A*a*b*c^2*d^4*tan(f*x + e) - 6*C*a*b*c^2*d^4*tan(f*x + e) - 3*B*b^2*c^2*d^4*tan
(f*x + e) - 2*A*a^2*c*d^5*tan(f*x + e) + 2*C*a^2*c*d^5*tan(f*x + e) + 4*B*a*b*c*d^5*tan(f*x + e) + 2*A*b^2*c*d
^5*tan(f*x + e) - B*a^2*d^6*tan(f*x + e) - 2*A*a*b*d^6*tan(f*x + e) + C*b^2*c^6 - C*a^2*c^4*d^2 - 2*B*a*b*c^4*
d^2 - A*b^2*c^4*d^2 + 3*C*b^2*c^4*d^2 + 2*B*a^2*c^3*d^3 + 4*A*a*b*c^3*d^3 - 4*C*a*b*c^3*d^3 - 2*B*b^2*c^3*d^3
- 3*A*a^2*c^2*d^4 + C*a^2*c^2*d^4 + 2*B*a*b*c^2*d^4 + A*b^2*c^2*d^4 - A*a^2*d^6)/((c^4*d^3 + 2*c^2*d^5 + d^7)*
(d*tan(f*x + e) + c)))/f

Mupad [B] (verification not implemented)

Time = 33.60 (sec) , antiderivative size = 3958, normalized size of antiderivative = 9.49 \[ \int \frac {(a+b \tan (e+f x))^2 \left (A+B \tan (e+f x)+C \tan ^2(e+f x)\right )}{(c+d \tan (e+f x))^2} \, dx=\text {Too large to display} \]

[In]

int(((a + b*tan(e + f*x))^2*(A + B*tan(e + f*x) + C*tan(e + f*x)^2))/(c + d*tan(e + f*x))^2,x)

[Out]

(log((2*C^2*b^4*c^5 - 2*C^2*a^2*b^2*c^5 + 4*C^2*b^4*c^3*d^2 - A*B*a^4*d^5 - 2*A*C*b^4*c^5 + B*C*a^4*d^5 + 2*A^
2*a*b^3*d^5 - 2*A^2*a^3*b*d^5 - A^2*a^4*c*d^4 + 2*B^2*a^3*b*d^5 - A^2*b^4*c*d^4 + B^2*a^4*c*d^4 + B^2*b^4*c*d^
4 - C^2*a^4*c*d^4 + C^2*b^4*c*d^4 - 4*C^2*a^2*b^2*c^3*d^2 + 5*A*B*a^2*b^2*d^5 + 2*A*C*a^2*b^2*c^5 + A*B*a^4*c^
2*d^3 + 3*A*B*b^4*c^2*d^3 - B*C*a^2*b^2*d^5 - 4*A*C*b^4*c^3*d^2 - B*C*a^4*c^2*d^3 - 3*B*C*b^4*c^2*d^3 + 2*B^2*
a*b^3*c^4*d - 2*C^2*a*b^3*c^4*d + 2*C^2*a^3*b*c^4*d - 2*A^2*a*b^3*c^2*d^3 + 6*A^2*a^2*b^2*c*d^4 + 2*A^2*a^3*b*
c^2*d^3 + 6*B^2*a*b^3*c^2*d^3 - 6*B^2*a^2*b^2*c*d^4 - 2*B^2*a^3*b*c^2*d^3 - 6*C^2*a*b^3*c^2*d^3 + 4*C^2*a^2*b^
2*c*d^4 + 6*C^2*a^3*b*c^2*d^3 - 2*A*C*a*b^3*d^5 + 2*A*C*a^3*b*d^5 - 4*B*C*a*b^3*c^5 + A*B*b^4*c^4*d + 2*A*C*a^
4*c*d^4 - B*C*b^4*c^4*d - 8*A*B*a*b^3*c*d^4 + 8*A*B*a^3*b*c*d^4 + 2*A*C*a*b^3*c^4*d - 2*A*C*a^3*b*c^4*d + 4*B*
C*a*b^3*c*d^4 - 8*B*C*a^3*b*c*d^4 - A*B*a^2*b^2*c^4*d + 8*A*C*a*b^3*c^2*d^3 - 10*A*C*a^2*b^2*c*d^4 - 8*A*C*a^3
*b*c^2*d^3 - 8*B*C*a*b^3*c^3*d^2 + 5*B*C*a^2*b^2*c^4*d - 8*A*B*a^2*b^2*c^2*d^3 + 4*A*C*a^2*b^2*c^3*d^2 + 16*B*
C*a^2*b^2*c^2*d^3)/(d^2*(c^2 + d^2)^2) + ((a*1i - b)^2*((A*b^2*d^2 - A*a^2*d^2 + C*a^2*d^2 - 8*C*b^2*c^2 - C*b
^2*d^2 + 2*B*a*b*d^2 + 4*B*b^2*c*d + 8*C*a*b*c*d)/d - (tan(e + f*x)*(3*B*a^2*d^5 - 5*B*b^2*d^5 - 4*C*b^2*c^5 +
6*A*a*b*d^5 - 10*C*a*b*d^5 + 4*A*a^2*c*d^4 - 4*A*b^2*c*d^4 + 2*B*b^2*c^4*d - 4*C*a^2*c*d^4 + 8*C*b^2*c*d^4 -
B*a^2*c^2*d^3 + B*b^2*c^2*d^3 - 8*B*a*b*c*d^4 + 4*C*a*b*c^4*d - 2*A*a*b*c^2*d^3 + 2*C*a*b*c^2*d^3))/(d^2*(c^2
+ d^2)) + (d*(a*1i - b)^2*(4*c*d - c^2*tan(e + f*x) + 3*d^2*tan(e + f*x))*(A + B*1i - C)*1i)/(c*1i - d)^2)*(A
+ B*1i - C)*1i)/(2*(c*1i - d)^2) + (tan(e + f*x)*(A^2*a^4*d^5 + A^2*b^4*d^5 + B^2*b^4*d^5 + C^2*a^4*d^5 + C^2*
b^4*d^5 - 2*A^2*a^2*b^2*d^5 + 3*B^2*a^2*b^2*d^5 + B^2*a^4*c^2*d^3 + 2*C^2*a^2*b^2*d^5 + 3*B^2*b^4*c^2*d^3 - 2*
A*C*a^4*d^5 - 2*A*C*b^4*d^5 - 2*B*C*b^4*c^5 - 4*C^2*a*b^3*c^5 + B^2*b^4*c^4*d + 4*A^2*a^2*b^2*c^2*d^3 - 4*B^2*
a^2*b^2*c^2*d^3 + 12*C^2*a^2*b^2*c^2*d^3 + 2*B*C*a^2*b^2*c^5 - 4*B*C*b^4*c^3*d^2 + 4*A^2*a*b^3*c*d^4 - 4*A^2*a
^3*b*c*d^4 - 4*B^2*a*b^3*c*d^4 + 4*B^2*a^3*b*c*d^4 - 4*C^2*a^3*b*c*d^4 - B^2*a^2*b^2*c^4*d - 8*C^2*a*b^3*c^3*d
^2 + 4*C^2*a^2*b^2*c^4*d + 2*A*B*a*b^3*d^5 - 4*A*B*a^3*b*d^5 + 4*A*C*a*b^3*c^5 - 2*A*B*a^4*c*d^4 - 2*A*B*b^4*c
*d^4 + 2*B*C*a^3*b*d^5 + 2*B*C*a^4*c*d^4 - 2*A*B*a*b^3*c^4*d - 4*A*C*a*b^3*c*d^4 + 8*A*C*a^3*b*c*d^4 + 4*B*C*a
*b^3*c^4*d - 2*B*C*a^3*b*c^4*d - 8*A*B*a*b^3*c^2*d^3 + 12*A*B*a^2*b^2*c*d^4 + 4*A*B*a^3*b*c^2*d^3 + 8*A*C*a*b^
3*c^3*d^2 - 4*A*C*a^2*b^2*c^4*d + 12*B*C*a*b^3*c^2*d^3 - 10*B*C*a^2*b^2*c*d^4 - 8*B*C*a^3*b*c^2*d^3 - 16*A*C*a
^2*b^2*c^2*d^3 + 4*B*C*a^2*b^2*c^3*d^2))/(d^2*(c^2 + d^2)^2))*(A*b^2 - A*a^2 - B*a^2*1i + B*b^2*1i + C*a^2 - C
*b^2 - A*a*b*2i + 2*B*a*b + C*a*b*2i))/(2*f*(2*c*d - c^2*1i + d^2*1i)) + (log((2*C^2*b^4*c^5 - 2*C^2*a^2*b^2*c
^5 + 4*C^2*b^4*c^3*d^2 - A*B*a^4*d^5 - 2*A*C*b^4*c^5 + B*C*a^4*d^5 + 2*A^2*a*b^3*d^5 - 2*A^2*a^3*b*d^5 - A^2*a
^4*c*d^4 + 2*B^2*a^3*b*d^5 - A^2*b^4*c*d^4 + B^2*a^4*c*d^4 + B^2*b^4*c*d^4 - C^2*a^4*c*d^4 + C^2*b^4*c*d^4 - 4
*C^2*a^2*b^2*c^3*d^2 + 5*A*B*a^2*b^2*d^5 + 2*A*C*a^2*b^2*c^5 + A*B*a^4*c^2*d^3 + 3*A*B*b^4*c^2*d^3 - B*C*a^2*b
^2*d^5 - 4*A*C*b^4*c^3*d^2 - B*C*a^4*c^2*d^3 - 3*B*C*b^4*c^2*d^3 + 2*B^2*a*b^3*c^4*d - 2*C^2*a*b^3*c^4*d + 2*C
^2*a^3*b*c^4*d - 2*A^2*a*b^3*c^2*d^3 + 6*A^2*a^2*b^2*c*d^4 + 2*A^2*a^3*b*c^2*d^3 + 6*B^2*a*b^3*c^2*d^3 - 6*B^2
*a^2*b^2*c*d^4 - 2*B^2*a^3*b*c^2*d^3 - 6*C^2*a*b^3*c^2*d^3 + 4*C^2*a^2*b^2*c*d^4 + 6*C^2*a^3*b*c^2*d^3 - 2*A*C
*a*b^3*d^5 + 2*A*C*a^3*b*d^5 - 4*B*C*a*b^3*c^5 + A*B*b^4*c^4*d + 2*A*C*a^4*c*d^4 - B*C*b^4*c^4*d - 8*A*B*a*b^3
*c*d^4 + 8*A*B*a^3*b*c*d^4 + 2*A*C*a*b^3*c^4*d - 2*A*C*a^3*b*c^4*d + 4*B*C*a*b^3*c*d^4 - 8*B*C*a^3*b*c*d^4 - A
*B*a^2*b^2*c^4*d + 8*A*C*a*b^3*c^2*d^3 - 10*A*C*a^2*b^2*c*d^4 - 8*A*C*a^3*b*c^2*d^3 - 8*B*C*a*b^3*c^3*d^2 + 5*
B*C*a^2*b^2*c^4*d - 8*A*B*a^2*b^2*c^2*d^3 + 4*A*C*a^2*b^2*c^3*d^2 + 16*B*C*a^2*b^2*c^2*d^3)/(d^2*(c^2 + d^2)^2
) + (tan(e + f*x)*(A^2*a^4*d^5 + A^2*b^4*d^5 + B^2*b^4*d^5 + C^2*a^4*d^5 + C^2*b^4*d^5 - 2*A^2*a^2*b^2*d^5 + 3
*B^2*a^2*b^2*d^5 + B^2*a^4*c^2*d^3 + 2*C^2*a^2*b^2*d^5 + 3*B^2*b^4*c^2*d^3 - 2*A*C*a^4*d^5 - 2*A*C*b^4*d^5 - 2
*B*C*b^4*c^5 - 4*C^2*a*b^3*c^5 + B^2*b^4*c^4*d + 4*A^2*a^2*b^2*c^2*d^3 - 4*B^2*a^2*b^2*c^2*d^3 + 12*C^2*a^2*b^
2*c^2*d^3 + 2*B*C*a^2*b^2*c^5 - 4*B*C*b^4*c^3*d^2 + 4*A^2*a*b^3*c*d^4 - 4*A^2*a^3*b*c*d^4 - 4*B^2*a*b^3*c*d^4
+ 4*B^2*a^3*b*c*d^4 - 4*C^2*a^3*b*c*d^4 - B^2*a^2*b^2*c^4*d - 8*C^2*a*b^3*c^3*d^2 + 4*C^2*a^2*b^2*c^4*d + 2*A*
B*a*b^3*d^5 - 4*A*B*a^3*b*d^5 + 4*A*C*a*b^3*c^5 - 2*A*B*a^4*c*d^4 - 2*A*B*b^4*c*d^4 + 2*B*C*a^3*b*d^5 + 2*B*C*
a^4*c*d^4 - 2*A*B*a*b^3*c^4*d - 4*A*C*a*b^3*c*d^4 + 8*A*C*a^3*b*c*d^4 + 4*B*C*a*b^3*c^4*d - 2*B*C*a^3*b*c^4*d
- 8*A*B*a*b^3*c^2*d^3 + 12*A*B*a^2*b^2*c*d^4 + 4*A*B*a^3*b*c^2*d^3 + 8*A*C*a*b^3*c^3*d^2 - 4*A*C*a^2*b^2*c^4*d
+ 12*B*C*a*b^3*c^2*d^3 - 10*B*C*a^2*b^2*c*d^4 - 8*B*C*a^3*b*c^2*d^3 - 16*A*C*a^2*b^2*c^2*d^3 + 4*B*C*a^2*b^2*
c^3*d^2))/(d^2*(c^2 + d^2)^2) + ((a*1i + b)^2*((tan(e + f*x)*(3*B*a^2*d^5 - 5*B*b^2*d^5 - 4*C*b^2*c^5 + 6*A*a*
b*d^5 - 10*C*a*b*d^5 + 4*A*a^2*c*d^4 - 4*A*b^2*c*d^4 + 2*B*b^2*c^4*d - 4*C*a^2*c*d^4 + 8*C*b^2*c*d^4 - B*a^2*c
^2*d^3 + B*b^2*c^2*d^3 - 8*B*a*b*c*d^4 + 4*C*a*b*c^4*d - 2*A*a*b*c^2*d^3 + 2*C*a*b*c^2*d^3))/(d^2*(c^2 + d^2))
- (A*b^2*d^2 - A*a^2*d^2 + C*a^2*d^2 - 8*C*b^2*c^2 - C*b^2*d^2 + 2*B*a*b*d^2 + 4*B*b^2*c*d + 8*C*a*b*c*d)/d +
(d*(a*1i + b)^2*(4*c*d - c^2*tan(e + f*x) + 3*d^2*tan(e + f*x))*(A*1i + B - C*1i))/(c*1i + d)^2)*(A*1i + B -
C*1i))/(2*(c*1i + d)^2))*(A*b^2*1i - A*a^2*1i - B*a^2 + B*b^2 + C*a^2*1i - C*b^2*1i - 2*A*a*b + B*a*b*2i + 2*C
*a*b))/(2*f*(c*d*2i - c^2 + d^2)) - (log(c + d*tan(e + f*x))*(d^3*(B*a^2*c^2 - 3*B*b^2*c^2 + 2*A*a*b*c^2 - 6*C
*a*b*c^2) - d^5*(B*a^2 + 2*A*a*b) - d*(B*b^2*c^4 + 2*C*a*b*c^4) + d^4*(2*A*b^2*c - 2*A*a^2*c + 2*C*a^2*c + 4*B
*a*b*c) + 2*C*b^2*c^5 + 4*C*b^2*c^3*d^2))/(f*(d^7 + 2*c^2*d^5 + c^4*d^3)) + (C*b^2*tan(e + f*x))/(d^2*f) - (A*
a^2*d^4 + C*b^2*c^4 - B*a^2*c*d^3 - B*b^2*c^3*d + A*b^2*c^2*d^2 + C*a^2*c^2*d^2 - 2*A*a*b*c*d^3 - 2*C*a*b*c^3*
d + 2*B*a*b*c^2*d^2)/(d*f*(c*d^2 + d^3*tan(e + f*x))*(c^2 + d^2))